clear;clc;clf;
syms f(x)
f=sin((1025/1000)*x)^2;
fplot(f,[-5,5])
grid on;

e=linspace(-5,5,100);  %划分区间
p=subs(f,x,e);  %区间中点的函数值
s_1=vpa(mean(p),20)   %计算平均值
 
d=linspace(-5,5,200);  %划分区间
p=subs(f,x,d);  %区间中点的函数值
s_2=vpa(mean(p),20)   %计算平均值

c=linspace(-5,5,1000);  %划分区间
p=subs(f,x,c);  %区间中点的函数值
s_3=vpa(mean(p),20)   %计算平均值

x_1=vpa(solve(f-s_3,x),20)
S=vpa(int(f,x,-5,5),20)